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The Graeffe Process as Applied to 
| Power Series 


Of the many methods which have been proposed for solving algebraic 
equations the most practical one, where complex roots are concerned, is the 
well known “root-squaring’”’ method usually referred to as the GRAEFFE! 
process. Nearly all accounts of the method show by examples how all the 
roots of an equation may be found by first finding all the absolute values 
of the roots, and later by various devices their arguments. The problem of 
finding the absolute values of the roots is comparatively simple and can 
be extended (with certain practical difficulties) to the case of the zeros of 
entire functions. This step was made early in the history of the problem 
by EULER? in generalizing Bernoulli's method to determine the first few 
zeros of the Bessel function Jo(z). More recent accounts of the problem of 
finding the absolute values of the zeros of entire functions have been given 
by Ostrowsk!! and Pérya.* The problem of finding by the Graeffe process 
the arguments of the zeros of entire functions (when these zeros are com- 
plex) apparently has not been attempted. At first this may seem a little 
strange since the Graeffe process owes its popularity to the fact that it is 
so successful in finding the complex roots of algebraic equations. The expla- 
nation is that the various devices, referred to above, for finding the argu- 
ments of the complex roots depend on an equation which has a finite 
number of roots. 

There is however one modification of the Graeffe process for poly- 
nomials, given more than a score of years ago by Bropetsky & SMEAL‘ 
and apparently overlooked by recent writers, in which a pair of complex 
roots may be found without having to find any other roots. The purpose 
of this note is to call attention to the fact that this process may be extended 
to entire functions and their complex zeros. 

As mentioned above, accounts of the Graeffe process are nearly always 
given in terms of examples. To be sure, the ‘general ideas’’ of the process 
are also set forth, but I have yet to see a writer admit that the Graeffe 
process fails utterly in such simple cases as the equation 


for the primitive fifth roots of unity. Although the Graeffe process is of 
considerable practical importance its theoretical justification, including the 
estimation of errors, is not at all easy. Even the lengthy memoir of Ostrowski, 
which is concerned largely with the absolute values of the roots, leaves 
something to be desired. If one is given only the coefficients of an equation 
with no advance information as to the disposition of its roots almost any- 
thing may happen if one applies the Graeffe process. Ostrowski shrewdly 
points out the distinction between the a priori and the a posteriori evalua- 
tion of the error, i.e. before and after the root-squaring process has been 
applied. It is clear that comparatively little of practical value may be said 
of the former. Even the number of steps necessary cannot be easily foreseen. 
The latter, unfortunately, depends on too many contingencies to afford a 
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378 THE GRAEFFE PROCESS AS APPLIED TO POWER SERIES 


simple treatment. In the brief discussion of the following method we must 
therefore confine ourselves to a few typical cases, leaving a more detailed 
and rigorous discussion to a later date and a more general method. 

Let 


(1) F(z) = dy + + + --- (a; real) 


be an entire function whose zeros are required. We may suppose that 
@> = 1. We proceed to calculate for m = 0,1, 2, --- two sequences of 
numbers 


defined as follows: For m = 0 
A, = a, B, = (r — 1)ay4 


so that Ap = 1, Bh = = 0. Form =1 


r—1 
= + 2D (—1)’a0e~ 
v=0 


r—1 
BLO = (—1)"(2r — 2v — 
v=0 
and for m > 1 


= + 25 (—1)"A 
v=0 


(3) ar 
= (—1)"A, BYR). 
v=0 
For all m: 
(4) Ao = 1, = 0. 


Theoretically, the zeros of f(z) are, in general, limits, as m— ©, of 
simple functions of A,“ and B,™. In what follows it will be convenient to 
let s denote 2”. Let the reciprocals of the zeros of (1) be denoted by 
a(R = 1,2, ---), the notation being such that 


(S) 2 |as| 2 |as| 2 ---. 


Finally let ¢ be an independent variable. Then the function defined recur- 
rently by 


a= al) = all) = — 


gm(t) = Vt) V2), 
has the expansion 


(6) en(t) = + + 0(2) 
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= (7) 
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THE GRAEFFE PROCESS AS APPLIED TO POWER SERIES 379 
where A,™, B,™ are given by (3). Moreover the zeros of this function 
are at 
(7) 


Identifying the elementary symmetric functions of (7) with the coefficients 
of (6) and recalling the definition of « and (4), we have, in the notation of 
symmetric functions 


(8) = > = > 


It is from these equations that we derive approximations to a,. Differ- 
ent cases arise, depending on which of the 2 signs in (5) may be replaced 
by =. The simplest case is that in which 


(9) > > |as| > ---. 


Although one might call this the ‘“‘general’’ case, a moment’s reflection 
shows that in this case all the zeros of f(z) are real. This is therefore a special 
case. By the usual argument, the symmetric functions (8) are, in this case, 
dominated by their first terms, provided s is so large that A,™ is sensibly 
the square of A,“-». Hence, approximately 


Ay™ — a? By™ = a; 


In general, equations of this type continue as long as > signs continue to 
hold in (9). A knowledge of the A’s alone suffices to determine |o,| approxi- 
mately by the usual relations 


| = | | 
(10) Jas] = | 
| | = | As | 


To determine the signs of these real roots it has been customary® to 
substitute +a, into the equation f(z) = 0, or to draw a graph of f(z). But 
knowledge of the B’s provides a rational method of determining the values 
of a, (with their signs). In fact, we have approximately 


a = Ay™ /By™ 


(11) = /By™ 
a3 = A;™ 


The agreement of the two values of |a,| obtained in (10) and (11) affords 
an excellent check on the work. If, however, the B’s are calculated only 
roughly, the signs of a, can be determined from the signs of B,™, since 
s — 1 is odd. 
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Suppose now that there are zeros of equal absolute value. The most 
important case of this kind is that in which the a’s occur in pairs of complex 
conjugate values, no two pairs having the same absolute value. This has 
been called the ‘general case.” To discuss it briefly we write 


a = a2 = ---, ani = = 
where the absolute values p; satisfy 
Pi > po > ps > 


For large s, (8) gives in this case, as dominant terms, 


Ayw™ = — cos sO; 
= 

= — cos 
= pip,” 


(12) ByY™ cos (s 1); 
= cos A 
= — cos (s — 1)62 + cos 6; cos 
= 2p;"*p2"*-! cos + cos 


This case is therefore recognized by the oscillatory character of ASP,; as 
functions of m, together with the fact that A{%? approaches the square of 
(k = 0,1, 2, ---). It is now clear how we may determine (p91, 
(p2, 92), -- + and verify our results. In fact 

= | | 1/26, = | /Ag™ | 1/20, oes 
(13) cos = /Ax™, 2.cos & = 


The values of (p2, 02), (ps, 3), thus found may be verified by 
using the expressions of A{”, A$”, --- and BY”, BY”, ---. If the B’s are 
calculated only roughly the values of 6;, 62, --- obtained from them will be 
sufficiently exact to determine which one of the s a priori values of 6 as 
defined by 


— ASP, 


= 
*  V[ASPAS? 


is the correct 6, to choose. The case of more than two zeros having equal 
absolute values breaks into subcases, each recognizable from the behavior 
of the A’s and each capable of a similar though more complicated treatment. 

Of course f(x) = 0 may have some real and some complex zeros. In 
practice this is especially true of polynomials but rather rare in the case of 
typical entire functions. In this case the above discussions are only slightly 
modified. Thus if the first three zeros ai, a2, a3 are real and all the other 
a’s are complex, equations (10) and (11) hold for a1, a2, and a3, while equa- 
tions (12) and (13) are modified by replacing Af”, BY” by Af?;/A$”, 
Bf?;/A$” in accordance with the fundamental equations (8). 

We consider now certain practical difficulties which arise when one 
extends the Graeffe process to entire functions. The most conspicuous of 
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THE GRAEFFE PROCESS AS APPLIED TO POWER SERIES 381 
these is due to the fact that the A,~™ and B,™ are each functions of the 
first 2”r coefficients of f(z) so that even for moderate values of r and m 
a large number of these coefficients must be known. Of course if f(z) is 
merely a polynomial of moderate degree the difficulty does not arise and 
most of the terms in the calculation of A,“ and B,™ vanish. To put it 
another way, the Graeffe process for a polynomial f of degree » makes 
considerable use of the fact that the coefficients of f, after the n-th, are 
exactly zero. This should warn the computer to be careful in truncating a 
power series if he expects the roots of the polynomial thus obtained to 
agree with the zeros of the series, especially for large zeros. In many typical 
cases the coefficients a, for m large may be closely approximated not by 
zero but by some simple asymptotic formula. This fact will reflect itself in 
the formulas (3) for A and B, thus enabling one to estimate the contribution 
from the far distant coefficients. In many cases also the function f(z) pos- 
sesses a simple asymptotic expansion for large z which enables one to 
compute (with bounded accuracy) a large zero of f(z). In these cases the 
Graeffe process may be used for finding only those zeros which are com- 
paratively near the origin. 

Another source of trouble is loss of accuracy in computing the A’s and 
B’s. This may be due to the accumulation of errors from “rounding off” 
the coefficients in the given equation. If enough accuracy can be main- 
tained until the process begins to converge, that is, until the expressions 
for A-™ and B,™ in (3) are dominated by their first few terms, then the 
crisis is passed and no great loss of accuracy occurs thereafter. This diffi- 
culty increases when the roots are nearly equal and is intolerable (as 
Ostrowski points out) when all the roots coincide, a case which is however 
of no practical importance. 

A more insidious source of trouble is a loss of accuracy due to cancella- 
tion of terms in computing A,“™ and B,“™. Here the extreme case is that 
in which f(z) = exp g(z), where g is an odd function. In this case one may 
show that A,“ = 0 for k > 0. In particular if f(z) = e** the B’s also vanish 
with the exception of B,“ = 2c. This puts an end to the Graeffe process. 
Of course in this case the function f(z) has no zeros anyway. The converse 
proposition is also true: If Ax = 0 for k > 0, and Bk = 0 for k > 1, 
then f(z) = e** and c = 3B,™. Hence in practice we need not expect an 
outright failure of the process. However, if the coefficients differ but slightly 
from those of an exponential function one may expect a considerable can- 
cellation of terms in computing A,™ and B,™, resulting in a devastating 
loss of accuracy. Hence extreme accuracy is required in the given coeffi- 
cients a,. If the a, are defined as rational numbers it may be advisable to 
use these exact values rather than any approximate decimal equivalents. 
In this way the A’s and B’s appear as exact rational numbers, no accuracy 
being lost. This device may have to be used for several stages until the 
process shows signs of converging. If the a, are given to within only a 
certain degree of accuracy and there is this tendency toward cancellation, 
very little can be done about it and some other method may have to be 
used. This difficulty is only slightly noticeable when f(z) is a polynomial of 
low degree. It is quite pronounced in the following example. 

Example. The following equation, proposed to me by F. Wolf, arises in 
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382 THE GRAEFFE PROCESS AS APPLIED TO POWER SERIRS 
connection with the flow of a hot gas through a pipe: 
J?(2) = Jo(z)J2(z) 


where J,(z) is the ordinary Bessel function of the first kind. It may be 
shown that 


2A(z) = 2[J2(z) — Jo(z)J2(z)] = 4 nJ,2(z) 
n=2 
from which it follows that except for the trivial double zero at z = 0, 


A(z) has no real or imaginary zeros. Replacing 2? by —4x and expanding 
Jo, J; and Jz in powers of x one finds that 


2m+2 x™ 
m 


If we wish to compute for example A{} and we carry the calculations of 
the terms to ten significant figures we obtain (in units of the 33rd decimal 
place) 


— ay = —28013 32771 
2a. 48133 66627 
ais = —30467 02393 
2as a= 14129 74452 
4754 97203 
2d¢ ae = 1144 16515 
—2as 192 74739 


2a, is = 22 06130 
—2a; 1 64295 
2a2 = 7456 
2a, aay => 185 
2ao ana = 2 


A? = = —404 


Hence only three ‘significant figures of A{? are obtained from 10-figure data. 
If however the a’s are taken at their rational values the rational value of 
A$) is found to be exactly 


—19471 33037 / (2-319. 58-75-114-13?-17-19-23) = 


—4,04020 96279 4... -10-. 
In this example the first 32 coefficients of f(x) were used to find 
A,;™ = 1.9683 98486 -10-” = —1.4572 73376 -10-" 
A, = 1.7561 21066 -10-*4 BL = —1.5624 51484 -10-% 
14 A, = —3.62350470 -10-% = 2p. cos 32 
(14) = 1.3060 43223 -10-* = 2% cos 316, 
A,® = —8.2527 36680 -10-7 B,® = —3.3783 23209 -10* 


We first find that 
= 0.18098 48935 
= —1.6102 54104. 
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These are approximations to p: and 2 cos 6; respectively. Substituting these 
into (14) we find by successive approximations that (14) has the simul- 
taneous solution 

pi = 0.1809 84894 

6, = 2.5066 8639. 


In this calculation it is not feasible to carry the main root-squaring process 
any further. One could have wished for more coefficients and higher powers, 
to separate more zeros. In this case the remaining zeros can be approxi- 
mated by an asymptotic formula with sufficient accuracy except in the 
case of a3 and a, which are still fairly close to the origin. The formula gives 
only 


p2 = 0.063, 4 = 2.70. 
Since the larger zeros are well determined and since we now know the first 
pair of complex zeros, it is a simple matter to compute the sums of the kth 
powers of all the zeros except a3; and a, for k = 7 and 8. From the above 
values of A; and B,® it follows that 
a3® + a? = 2° cos 86. = —5.4461-10-” 
as? + ay? = 2p? cos 70. = 9.1125-10-° 


The approximate values of p2 and 62 given above make it easy to determine 
by successive approximations the solution 


p2 = 0.064339, 6 = 2.70071 


of the above system. 


The above example is given in such detail in order to indicate what 
steps may be necessary in dealing with intractable entire functions. When 
one is deprived of the use of extremely high powers of the zeros (custom- 
arily used in treating polynomials) one is obliged to obtain more information 
than usual from the lower powers. 


D. H. L. 


1For a typical account of this method, together with some historical remarks, see 
E. T. WHITTAKER & G. Rosinson, The Calculus of Observations, third ed., London, 1940, 
p. 106-120. References to other accounts are given in C. A. Hutcurinson, “On Graeffe’s 
method for the numerical solution of algebraic equations,” Am. Math. Mo., v. 42, 1935, 
p. 149-161, and L. L. Cronvicn, “‘On the Graeffe method of solution of equations,’’ Am. 
Math. Mo., v. 46, 1939, p. 185-190. An early exponent of this method, who did much to 
make it well known, was the astronomer J. F. ENcKE (J. f. d. reine u. angew. Math., v. 22, 
1841, p. 193f). He attributed the method to Graeffe. He was also responsible for an utterly 
useless scheme of changing the signs of all the roots. Many writers on the subject still use 
this confusing notion of the ‘Encke roots” of an equation. Graeffe’s method of solutions 
of numerical equations was really due to DANDELIN (1794-1847), Mem. Acad. Brussels, 
v. 3, 1826, p. 46, or to WARING (1734-1798), Meditationes Analyticae, 1776, p. 311. It was 
published by C. H. Grarrre (1799-1873) in a prize paper, Die Auflésung der hiheren 
numerischen Gleichungen, Ziirich, 1837. Other accounts of the history of the method will 
be found in A. M. OstrowskI, “‘Recherches sur la méthode de Graeffe et les zéros des poly- 
nomes et des séries de Laurent,’’ Acta Math., v. 72, 1940, p. 99-155, and E. CarvALLo, 
Méthode pratique pour la Résolution numérique complete des équations algébriques ou transcen- 
dantes, Paris diss., 1890. 

2L. Eucer, Akad. Nauk, S.S.S.R., Leningrad, Acta, 1781, part 1, 1784, p. 170f. 

3G. Pérya, “Uber das Graeffesche Verfahren,” Z. Math. Phys., v. 63, 1915, p. 275-290. 

4S. Bropetsky & G. SMEAL, ‘On Graeffe’s method for complex roots of algebraic 
equations,’’ Camb. Phil. So., Proc., v. 22, 1924, p. 83f. 

5 A. C. AITKEN has proposed a root-cubing procedure in which real roots are deter- 
mined together with their signs; see Math. Gazette, v. 15, 1931, p. 490-491. 
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184{A, B, C, D, E, M].—Epwarp ScHWITZER ALLEN (1887— ), Six-Place 
Tables. A Selection of Tables of Squares, Cubes, Square Roots, Cube Roots, 
Fifth Roots and Powers, Circumferences and Areas of Circles, Common 
Logarithms of Numbers and of the Trigonometric Functions, the Natural 
Trigonometric Functions, Natural Logarithms, Exponential and Hyper- 
bolic Functions, and Integrals. Sixth ed., New York, McGraw-Hill, 1941. 
xxiii, 181 p. 10.7 X 17.7 cm. Leatherette $1.50. Compare RMT 49. 


Through revisions, including enlargements, and the addition of marginal index and 
flexible covers, in new editions (1925, 1929, 1931, 1935, 1941) since the original 1922 pocket 
volume, this work has attained to its present very useful form. At many places of the com- 
pilation, the thorough competence of the author is in evidence. In this latest edition three 
new tables are added, namely: (a) T. VIII, a four-place table of sine, cosine, tangent, 
cotangent, and log sin, log cos, log tan, for radian argument [0(.01)3.20; 4D], included at 
the suggestion of electrical engineers; (b) T. [X, log m!, = [1(1)100; 6D]; and (c) a four- 
place table of complete elliptic integrals K, E, for sin k = 0(1°)65°(0°.5)80°(0°.2)89°(0°.1) 
89°.9 added to T. XIII, containing 167 integrals. 

T. VI is a six-place table of all six of the natural trigonometric functions for each min- 
ute of the quadrant; T. XI, for e*, e~*, sinh x, cosh x, tanh x, x = [0(.01)1(.05)2.4(.1)5(1)10; 


mostly 4D], and T(z), x = (0(.01).99; 4D], and x = [0(.01)2.99; 4D]. Thus 


not all tables are six-place; T. III of circumferences and areas of circles with diameters 
from 1/64 to 100 is mostly 6S. The last Table, XIV, is devoted to mathematical constants. 
In the xxiii pages of the Introduction the theory of logarithms and explanation of tables 
are given, with exercises for solution. 

A. 


185[A, B, D, E, G, H, I].—Fritz Empe, Tables of Elementary Functions . . . 
With 83 figures. Tafeln elementarer Funktionen . . . mit 83 Textfiguren. 
Leipzig and Berlin, Teubner, 1940. Two title-pages. On the back of the 
second title-page is, in part, Copyright vested in the Alien Property Cus- 
todian, 1944, pursuant to law. Published and distributed in the Public 
Interest by authority of the Alien Property Custodian under License No. 
A-709. Lithoprinted by Edwards Brothers, Inc., Ann Arbor, Michigan, 
U. S. A., 1945. xii, 181 p. 17.7 XK 26.1 cm. $3.20. 


There was incidental reference to this work in MTAC, p. 108, Oct. 1943 (when no 
copy of the work was available in this country), in the course of reviews of the six editions 
or reprints then available of JAHNKE & EmpbeE, Tables of Functions. The first 78 p., devoted 
to Elementary Functions in the third edition of 1938, did not appear in the fourth or fifth 
editions but were reprinted in the sixth. The work under review was intended to be a revi- 
sion, and in various directions greatly to elaborate the 1938 material. Larger type is used 
and the type-page is about 25% larger than that of the Tables of Functions of JAHNKE & 
Empe. Instead of X sections in the original 78 pages there are now XVII sections in the 
172 pages. We shall now briefly indicate most of what is new. Section I, Powers, p. 1-19 
(old, p. 1-11); the new tables are of squares and reciprocals, tables of cubes and reciprocals 
of the squares, powers of 2, 3, 5, 7, 11, x4, (10x)#, (x/10)"8, x1/8, (10x)*/3, Section II, Factor 
Table, Factorials, etc., p. 20-25 (new): Factor table of 49-9997; n!, 1/n!, n!!, 1/n!!, 
n = 1(1)20; multiples of loge = M, x/M, 2x/x, xx/2. Section III, Auxiliary tables for 
computation with complex numbers, p. 26-33 (old, p. 12-19), slight rearrangements. 
Section IV, Quadratic equations, p. 34-37 (new), graphs and tables for the solution of 


384 


Ag 4 
for | 
ang) 
also 
amt 
int 
pow 
and 
tan 
and 
sinh 
and 
witl 
nun 
not 
is cl 
exce 
for 
ada 
pre] 
Em 
be pro 
mal 
stuc 
18¢ 
: Cor 
nH, 
Def 


RECENT MATHEMATICAL TABLES 385 


x? + px + ¢ = 0. Section V, Cubic equations, p. 38-47 (old, p. 20-30), much the same. 
Section VI, Equations of the fourth degree, p. 47-72 (new), elaborate graphs, and tables 
for finding real and compiex roots of biquadratic equations. Section VII, Circular functions, 
angles in degrees, p. 73-79 (new), tables of sin x, cos x, tan x, 2 sin x, V2 cos x, 1 — cos x; 
also of integrals; graphs; etc. Section VIII, Angles in quadrants, p. 80-89 (new), tables of 
amp. sin x, sinh cos x, cosh e**, tan x, tanh coth }xx, V2 sin x, 
V2 cos x, 1 — cosx, cosh }rx — 1; also graphs. Section IX, Hyperbolic functions, angles 
in radians, p. 90-115 (old, p. 52-60), various formulae, functions of a multiple sector, 
powers, derivatives, integrals, approximate values for small sectors, gudermannian; graphs 
and altitude chart; numerous tables include those of sin x, cos x, tanx, cos x, 
tan™ x, sinh x, cosh x, tanh x, sinh x, cosh x, e*, e*, 1 — cos x, cosh x — 1; tables of Inx 
and log x. Section X, Special functions, p. 116-130 (old, p. 36-39, 44-51), tables of e~*’, 
Planck’s radiation function, Langevin’s function cothx—1/x, ptanp, — cot p/p, 
— tan p/p, p in right angles; numerous graphs including Hayashi’s sinh x sin x, cosh x cos x, 
sinh x cos x, cosh x sin x. Section XI, Transcendental equations, p. 130-131 (old, p. 30-31). 
Section XII, Circular and hyperbolic functions of a complex variable, p. 132-145 (old, p. 
60-76), some revisions. Section XIII, The function [tan (i7r)] /i#r, p. 146-155 (new), graphs 
and table of the function. Section XIV, Chebyshev’s polynomials, p. 156-157 (new), graphs. 
Section XV, The function z‘, p. 158 (new), graph. Section XVI, Approximate calculations 
with polynomials, p. 159-162 (new), formulae, tables. Section XVII, Some remarks on 
numerical calculations, p. 163-172 (new). The tables on old pages 32-35 and 40-44 do 
not appear in the newer work. 

On pages 173-176 is a list of nearly 100 “‘Useful books for the computer.” 

In various tables the interval in the argument is either 1 or 2 or 5 units in the last 
figure. The differences printed are calculated throughout for the interval 1. “‘The interval 
is chosen so great that the maximal error which might arise by linear interpolation does not 
exceed .5 (10-*) of the function value. Further details on this point are to be found in the 
last chapter. This increase of the interval reduces the tables considerably, so that the values 
for many functions can be given in small space. The tables are not overstretched, but 
adapted to the demand of each function.” 

The work is in English and German throughout, the English translations having been 
prepared by G. W. O. Howe of the University of Glasgow. 

Emde felt that his “‘book will be useful to many who never refer to the 1938 Jahnke- 
Emde. The latter also include most of the students of Technical Colleges. Our book may 
prove a help in their mathematical instruction. The numerous figures to scale, of which 
many show the same function in various scales, will make the functions clearer. An attentive 
study of the variation of the interval in the several functions will lead to greater insight.” 


186[A, D].—H. M. Terri, & M. “Tables of numbers 
related to the tangent coefficients,’ Franklin Institute, J., v. 239, Jan. 
1945, p. 66-67. 16.5 K 24.1 cm. 

If tang = Tie + 
= 229(2% — 1)B,/2n. 


Correcting the erroneous relation in the text, the number H, should have been defined by 
nH,, = 2"(2" — 1)Bn. Hence 


HA, = 
Defining the number K, = 2(2% — 1)B,, we have K, = nT,/2"~, and 
nH, = 2°-1K,. 
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The chief tables before us give the values of H, and K,, (each an integer) for m = 1(1)19. 
Since there is a table of 7, for m = 1(1)30, in Peters and Stein’s Anhang, p. 88, of Peters’ 
Zehnstellige Logarithmentafel, v. 1, Berlin, 1922, the computation of the present tables 
for » = 1(1)19 was a rather trivial matter. There are also two small tables of H? and Kz, 
for p = 1(1)4, m = 2(1)5; these numbers are defined by the recurrence formulae 


— (p + 1)(2p + 1)H2** = 

where Hf = K = 1. 

In his Théorie des Nombres, Paris, 1891, v. 1, p. 251, Lucas refers to Gen = (—1)*K, as 

Genocchi numbers (A. Genoccat, Annali d. Sci. Matem., v. 3, 1852, p. 395 f.), though 

considered by EuLEerR (Institutiones Calculi Differentialis, v. 2, 1755, p. 486), and Lucas 

gives the first six numbers of the Terrill table. A ms. table of K, for n=110, by D. H. L., 

is in the Library of the American Mathematical Society. 
R.C. A. 


187[B, D, E, G, H, L].—P. R. E. Jannxe & F. Empe, Tables of Functions 
with Formulae and Curves |Funktionentafeln mit Formeln und Kurven . . .] 
Published and distributed in the public interest by authority of the Alien 
Property Custodian under license No. A-57, Fourth revised edition, New 
York, Dover Publications, 1780 Broadway, 1945. [xv], 306, 76 p. 15.4 
X 29.8 cm. $3.75. 


In MTAC, p. 106-109, 171 and 204 (corrigenda), 293-294, some information has been 
given regarding the five previous editions or reprints in 1909, 1933, 1938, 1941, and 1943. 
Now we have the sixth, with the correction of 396 errors in the fifth edition (or reprint), 
and the addition of a Supplementary Bibliography containing 43 titles. Material for correc- 
tion of the errors was mainly furnished, either directly or indirectly, by R. C. A., L. J. C., 
NYMTP, and J. C. P. Mitxer. In many cases these errors, on 56 pages, are only unit 
errors in the last decimal place. A list of the pages where the only known errors in the 1943 
volume occur, followed by their number on each page, may be found to be of service. They 
are as follows: 

9 (1), 10 (2), 34 (8), 35 (3), 39 (1), 59 (2), 62 (1), 65 (1), 66 (1), 69 (3), 70 (1), 72 (3), 
74 (1), 79 (1), 85 (3), 94 (1), 102 (2), 103 (1), 105 (1), 117 (1), 124 (13), 125 (40), 128 (2), 
164 (32), 166 (1), 168 (27), 177 (1), 183 (1), 204 (2), 205 (1), 206 (4), 212 (1), 216 (2), 
218 (1), 220 (1), 228 (15), 229 (6), 234 (8), 235 (63), 236 (35), 237 (5), 238 (42), 239 (10), 
246 (1), 247 (2), 263 (3), 266 (1), 300 (3), 301 (1), [2] (2), [4] (4), [5] (8), [6] (5), [7] (12), 
[9] (1), [56] (2). The notation [56], for example, refers to p. 56 in the second section, of 
76 p., in the volume. There are also changes 298(6), 299(2). 

In earlier pages of MTAC a number of these errors have been already listed. References 
to these, a discussion of other errors of more than a unit in the final figure, as well as a list 
of errors in the 1909 and 1933 editions, are given elsewhere in this issue, MTE 61, p. 391 f. 
Thus information here placed at the disposal of owners of all editions will enable them 
to increase the value of their editions as tools, while purchasers of the volume under re- 
view will have a greatly improved edition of an extraordinarily useful work. 

A. 


188[C, D].—Five-Figure Logarithm Tables containing Logarithms of Num- 
bers and Logarithms of Trigonometrical Functions with Arguments in 
Degrees and Decimals. Published for the Ministry oF Suppty by His 
Majesty’s Stationery Office, London, 1944. iv, 74, XXII, 29-120 p. 
15 X 24 cm. 7s 6d. 


“This collection of 5-figure tables of logarithms and logarithmic trigonometrical func- 
tions has been compiled as a war-time measure, primarily to meet the urgent requirements 
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of the optical industry. Expediency can be the only justification for the juxtaposition of 
differing page sizes, of differing styles of types and of differing forms of printing.’’ So begins 
the Preface by H. J. Goucu, director-general of Scientific Research and Development, 
Ministry of Supply, writing in September 1944. 

The first section of the volume is a reproduction of p. 2-73 of E. CHapre., Five-Figure 
Mathematical Tables . . . , Edinburgh and London, W. & R. Chambers, 1915. Through 
an arrangement with the publishers the original stereotype plates were used. These 
reproduced pages give logarithms of numbers 10 000 to 40 000, and 4 000 to 10 000, with 
differences and proportional parts. See MTAC, p. 195. The next two parts are photolitho- 
graphic reproductions of German works, under license from the Custodian of Enemy 
Property. The first of these has the title-page: Die Logarithmen der Sinus und Tangenten 
fiir 0° bis 5° und der Cosinus und Cotangenten fiir 85° bis 90° von Tausendstel zu Tausendstel 
Grad als Ergiénzung zu C Bremiker’s 5 stelligen Logarithmentafeln, herausgegeben von M. 
von Rour . . . Berlin, Weidmannsche Buchhandlung, 1900, and the table occupies pages 
I-XX. The final part of the volume has the title page: Dr. C. Bremikers Logarithmisch- 
Trigonometrische Tafeln mit fiinf Decimalstellen, besorgt von . . . A. Katitus. Sechzehnte 
Stereotyp-Auflage, Berlin, Weidmannsche Buchhandlung, 1925. The pages 30-119 contain 
five-figure logarithms of sines, tangents, cotangents and cosines, for the angles 0(0°.01)45°. 
On the final page are values of constants. 

The original edition of these tables of Cart BREMIKER (1804-1877) appeared in 1872. 
Beginning with the third edition in 1880, A. Kallius was the editor. There was a seven- 
teenth edition in 1937. The third unchanged edition, of the supplement by Louis Otto 
Moritz von Rohr (1868—_ ), appeared in 1933. 

Explanations of any kind in connection with the tables were naturally regarded as 
unnecessary. Some American libraries will doubtless be glad of the opportunity offered for 
purchasing this volume so as to get a copy of Rohr’s table. 

Workers in the optical industry in this country, where computing machines are so 
freely used, are not likely to find any use for the volume under review. Rather would they 
turn to such a volume as the seven-place table of natural trigonometric functions for every 
thousandth of a degree, prepared by Peters in 1918 for the Goerz optical establishment, 
now readily available in an American edition. See MTAC, p. 12-13. 

R.C.A, 


189[D].—HERMANN BRANDENBURG, Sechsstellige trigonometrische Tafel alter 
Kreisteilung fiir Berechnungen mit der Rechenmaschine, enthaltend die 
unmittelbaren oder natiirlichen Werte der vier Winkellinien-Verhdltnisse: 
Sinus, Tangens, Cotangens und Cosinus des in 90° und 60’ geteilten 
Eiinheits- Viertelkreises in Unterschieden von 10 zu 10 Sekunden nebst zwei 
Vortafeln mit Einzelsekundenwerten der Cotangente und erweiterten Zehn- 
sekundenwerte des Sinus und der Tangente kleiner Winkel. . . . Das 
Vorwort in deutscher, englischer, franzdsischer, spanischer, hollindischer, 
russischer, und japanischer Sprache, Leipzig, Lorenz, 1932. Back of the 
title page is, in part, Copyright vested in the Alien Property Custodian, 
1944, pursuant to law. Published and distributed in the Public Interest by 
authority of the Alien Property Custodian under License No. A.407. Litho- 
printed by Edwards Brothers, Inc., Ann Arbor, Michigan, U. S. A., 1945. 
xxii, 304 p. 20.2 X 28 cm. $5.00. 


In 1923 Brandenburg published seven-place trigonometric tables containing natural 
values of the Sine, Tangent, Cotangent, and Cosine at interval 10”, with differences and 
proportional parts. There were also an introductory table of Cotangent up to 6°, and 
of Tangent from 84° to 90°, both at interval 1”. The construction of this work was based 
on the 10-place trigonometric table of RuatTicus, 1596, 1607, and the 15-place trigonometric 
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table of Prriscus, 1613 (see MTAC, p. 8-9, and A. DEMorcan, English Cyclopaedia, ed. 
C. Knight, Arts and Sciences, v. 7, London, 1861, article ‘““Table’”’). But this edition of 
Brandenburg’s work contained over 400 errors. All of these were eliminated in the second 
edition of 1931, which is “practically free from error.”” This work was prepared chiefly for 
computers working with calculating machines. 

From the material of these editions Brandenburg prepared six-place tables by rounding 
off the last decimal, and thus appeared the original of the work now under review. Vor- 
tafel 1 is of Cotangent to 3° (or of Tangent 87° to 90°) at interval 1”, to 7S; and Vortafel 
2 is of Sine and Tangent to 1°, at interval 10’, mostly to 7 or 8D. The last five pages con- 
tain miscellaneous numerical data and tables. Nine errata discovered in this volume by 
L. J. C. were published in MTAC, p. 162. Six of these, and two other corrections on p. 64, 
are listed on the back of the dedication page (p. iv). The two corrections not previously 
listed in MTAC are as follows: 


Diff. following cot 5° 42’ 30’, for 4903, read 4898 - 
P.P. for 9/10 of 490, for 441.0, read 441,0. 


The volume is well reproduced, and bound in light-blue buckram, with title in black ink 
along the back. 

Brandenburg well remarks that a ‘‘6-figure table has one advantage that must not be 
underestimated, namely that it suffices for the greater part of the calculations encountered 
in the solution of triangles, and because of its convenience in manipulation is particularly 
advantageous when five or fewer figures are sufficient, in which case the unused figures are 
rejected. Further the owner of a 7-figure table generally equips himself with a 6-figure 
table also, whereas the converse is rarely the case. The 6-figure table is undoubtedly espe- 
cially beloved in the computing profession, and can be pronounced the favorite among 
trigonometrical tables.” 

R. C. A, 


190[F].—L. ALaoctu & P. Erpdés, “On highly composite and similar 
numbers,” Amer. Math. So., Trans., v. 56, Nov. 1944, p. 448-469. 
19.5 X 27.5 cm. 


This paper contains two tables (p. 467-469) having to do with so-called abundant 
numbers. If the sum o(m) of the divisors of m exceeds 2n, then n is called abundant. In case 
o(n)>o(m) for all m<n, then n is called highly abundant. In case o(n)/n>o(m)/m for all 
m<xn, then n is called superabundant. 

The first table lists all highly abundant numbers < 10080. This table errs by including 
the non-highly abundant number 2280 and by excluding the highly abundant numbers 
1920 and 2160. The number of highly abundant numbers < 10080 is thus 84, not 83. 
The table also contains the minor errata: 


n = 1080 for 23.3.5, read 2.33.5 
n = 4620 for o(m) = 16120, read 16128. 


This table was based on Glaisher’s tables! of o(m) for n < 104. 

The second table lists all superabundant numbers < 10" in their factored form to- 
gether with values to 3D of o(m)/n in each case. These numbers are similar to Ramanujan’s? 
highly composite numbers which have more divisors than any smaller number. 

D. 


1J. W. L. GuatsHErR, Table I in BAASMTC, Mathematical Tables, volume VIII, Num- 
ber-Divisor Tables, Cambridge, 1940. 

2S. RAMANUJAN: “Highly composite numbers,”’ London Math. So., Proc., s. 2, v. 14, 
1915, p. 347-409; Collected Papers, Cambridge, 1927, p. 78-128. 
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191[(L].—R. P. BELL, “‘Eigen-values and eigen-functions for the operator 
@*/dx* — |x|,” Phil. Mag., s. 7, v. 35, Sept. 1944, p. 584-585. 16.9 x 
25.3 cm. 


T. I gives values of Jejs(z) and J_2(z), for z = [0(.02)1; 4D]. These were calculated 
from the following series: 


Jen(z) = .69783277(1 — .1500z* + .00703z* — .000162* + ---), and 
J-2ix(2) = .59255277/8(1 — .75002* + .070312* — .002512* + .0000472 — ---). 


Interpolation in T. I gave the first zero of J23(z) — J_2(z) as z = .6856, not quite agreeing 
with the value given by H. M. Macponatp, R. So. London, Proc., v. 90A, 1914, p. 54, 


where we find 16854, The eigen-values of (= |x| +r) ¥ = Oare given by 


— = 0, 2 even, 
— = 0, odd. 


The first six values of X,, are given, m = [0(1)5; 3D], in T. II. 


In T. III are the eigen-values A, of - | xt] +r)y=0 for =1,2, 4, 
nm = [0(1)5; 2 or 3D]. 
The error in Watson’s Bessel Functions, p. 751, noted by J. C. P. Miter, MTAC, 
p. 367, was already published by Bell, on p. 584. 
© a, 


192([L].—R. CLark Jones, “On the theory of the directional patterns of 
continuous source distributions on a plane surface,” Acoustical So. 
Amer., J., v. 16, Jan. 1945, p. 158-171. 19.8 X 26.6 cm. 


Computations are made for distributions of finite extent with a distribution function 
of type 
A(u) = (1 — u™)* 


with 10 different choices of (m, m) for both the circular and the linear cases. These choices 
are: 


(1, —1), (1,0), (1,1), (1,2), (1,3), (1,4), (2,1), (2, 2), (3, 1), (4, 1). 


There are the following 21 tables: T. I.—Jo(a) and J3(a) to 5S, 10log J>*(a) to 4S, 
for a = 0(.4)20. Also first 12 values of a for maxima of J3?(a) (zeros of Ji(a)) to 3D. 
T. Il.—Qy = 2Ji(a)/a and Q3, to 6S, 10log Qy* mostly 4S, a = 0(.1)20. Also first 12 
values of a for maxima of (zeros of J2(a)) to 3D. T. II].—Qw = 8J2(a)/a? and 
to 5S, and 10log Qy” to 4S, a = 0(.4)20. Also first 11 values of a for maxima of Q3, (zeros 
of Ji(a)) to 3D. T. IV.—Qy = 48Js(a)/a* and Q3, to 5S, and 10log Qy? to 2D, a = 0(.4)20. 
Also first 11 values of a for maxima of Q3, (zeros of Ji(a)), to 3D. T. V.—Qw = 384Js(a) /axt 
and Q%, to 5S, and 10log Qy? to 2D, a = 0(.4)20. Also first 10 values of a for maxima of 
Q3, (zeros of Js(a)), to 3D. T. VI.—Qwn = 3840Js(a)/a® and Q3, to 5S, and 10log Qy? to 4S, 
a = 0(.4)20. Also first 10 values of a for maxima of Q3, (zeros of Js(a)). T. VIIL.— 
Qv = 12J2(a)/a? — 24J3(a)/a* and Q3, to 5S, and 10log Qy” to 4S, a = 0(.4)20. T. VIII.— 
Qu = 120Js(a)/a* — 720Js(a)/a* + 1440J;(a)/a> and Q% to SS, and 10log to 4S, 
a = 0(.4)20. T. IX.—Qy = 16J2(a)/a* — 64J3(a)/a® + 128J,(a)/a*, and Q3, to 4S, and 
10log Qy* to 4S, a = 0(.4)20. T. X.—Qw = 20J2(a)/a® — 120J3(a)/a* + 
— 960J;(a)/o® and Qj, to 4S, and 10log Qy* to 4S. T. XI.—Qy = sin 6/8 to 6S, and Q}, 
to 5S, and 10log Qy? to 5S, 8 = 0(.4)20. Also first 5 values of 8, to 2D, for maxima of Q%,. 
T. XII.—Qw = 3(x/28)*Ja(8)/8 to 6S, Q2 to SS, and 10log Qx? to 5S, 8 = 0(.4)20. Also 
first 5 values of 8, 2-3D, for maxima of Q. T. XIII1.—Qw = 15(#/28)*Js2(8)/8* to 6S, 
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and Q3, to 5S, and 10log Qy” to 4S, 8 = 0(.4)20. Also first 4 values of 8, to 2D, for maxima 
of Q3,. T. XIV.—Qw = to 6S, to 5S, 10log Ox? to 4S, = 0(.4)20. 
Also first 4 values of 8, to 2D, for maxima of Q3,. T. XV.—Qw = 945(2/28)*J9/2(8)/B* to 6S, 
Q3, to 5S, 10log to 4S, = 0(.4)20. T. XVI.—Qw = — 
and Q3 to 5S, 10logQy? to 4S, @ = 0(.4)20. T. XVII.—Qw = (w/28)8[45Js(8)/B" 
— 270Jan(8)/6* + and to SS, and 10log to 4S, = 0(.4)20. T, 
XVIII.—Qw = — + 56Jr2(8)/6*], and Qj to 5S, 10log Qy* 
to 4S, 6 = 0(.4)20. T. XIX.—Qy = — + 216J72(8) /B* 
— 432Jn(6)/84], and to 5S, and 10log Qy? to 4S, = 0(.4)20. T. XXI.—Qy = 
cos — (28/x)*] and Q3, to SS, 10log Qx? to 4S, 8 = 0(.4)20. T. XXII.—Qy = sin 
— (8/x)?)] and Q3, to 5S, 10log Oy? to 4S, 8 = 0(.4)20. 
R. C. A. & H. B. 


193[L].— NATIONAL DEFENSE RESEARCH COMMITTEE, Division 6, Scatter- 
ing and Radiation from Circular Cylinders and Spheres, prepared by 
A. N. Lowan of the N. Y. M. T. P., & by P. M. Morse, H. FEsHBAcH, 
M. Lax, of the M. I. T. Underwater Sound Laboratory. Report Febru- 
ary 1945. iii + 124 leaves, 21.5 X 28 cm. Printed on only one side of 
each leaf, by the photo-offset process from manuscript. These tables are 
available only to certain Government agencies and activities. 


In problems of wave motion external to spheres and cylinders the Bessel functions of 
the first and second kinds are of importance. The tables 1-5, 8-12 in this report are all 
for x = 0(.1)10 with 2D, but the number of significant figures generally varies from 3 to 5. 
In tables 6, 7, 13 and 14, 6S are given and the ranges are x = 0(.1)10, m = 0(1)20. The 
same range for m is used for all the other tables except 10 and 11 where m = 0(1)15. The 
quantities tabulated are: 


Cylindrical Functions Spherical Functions 
1. an(x) = [— 8. an(x) = tan [— xjn’(x)/jn(x)] 
2. Bn(x) = tan (— xNn'(x)/Nn(x)] 9. Bn(x) = tan [— xmn’(x) /nn(x)] 
3. 5n(x) = tan [— Jn(x)/Na(x)] 10. 5,(x) = tan™ [— 7,(x)/mn(x)] 
4. 8,'(x) = [— Jn’ (x)/Nn'(x)] 11. 5,’(x) = [— jn’(x)/tn’(x)] 
5. n(x) = tan [tan 5, sec an Cos Bn] 12. ya(x) = tan [tan 6, sec an cos Bn] 
6. Cr(x) = |In(x) + iNn(x) | 13. Da(x) = |jn(x) + ma(x)| 
7. = |In'(x) + iNn’(x)| 14. = |jn'(x) + 


The function N,,(x) is the Bessel function of the second kind denoted by Y,(x) in MTAC. 
The functions j(z), (z) are defined by the equations 


jn(2) = = (4 


A useful list of properties of the functions is given. In the table of contents the expressions 
for the spherical functions 6,(x) and 6,’(x) in 10-11 are given incorrectly. 


194{L].—[a. M. Seresrtisxii, ‘“Obtekanie tel vrashcheniia” [Flow past 
bodies of revolution], Akad. Nauk, SSSR, Otdelenie tekhnichnikh Nauk, 
Institut Mekhaniki, Prikladnata Matematika i Mekhanika, Applied 
Mathematics and Mechanics, v. 8, June 1944, p. 103, 104. 16.7 & 25.5 cm. 


There is a table of Legendre functions Pa(u), dPn(u)/dy, for m = 1(1)5, w = [— 1(.05) 
—.9(.1)+.9(.05)1; 4D], where » = cos 0; @ is given to the nearest 1’ for each value of yu. 

There are also graphs of Q,(A) = Qn(cosh y), and of 5,(y), for y = .1(.05).35, m = 1(1)5, 
where dQ,,(A)/dA = — 


H. B. 
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Reference has been made to Errata under RMT187 (Jahnke & Emde), 
189 (Brandenburg), 190 (Alaoglu & Erdés), 191 (Watson), 193 (N. D. R. C.); 
N30 (Girard), 36 (Jacobi). 


61. E. JAHNKE & F. Empe, Tables of Functions. Editions as in RMT 113, 
p. 106f. The errors here listed are almost wholly supplementary to those 
discussed p. 108-109, 136-137, 161, 196, 198, 204 (corrigenda, p. 109, 
161), 286, 287, 294, 308 (corrigenda, p. 137), 325. 


The popularity on both sides of the Atlantic of the Dover Publications reprint (1943) 
of these well-known tables makes a list of their errors timely. The following lists do not 
pretend to be complete, but as they put on record the results of systematic examination of 
some of the tables, they may well relieve anxieties about the correctness of certain tables. 
I wish to place on record the great share that J. C. P. MILLER has had in these examinations, 
and in perfecting this list. While references have been previously made in MTAC to con- 
tributions of NYMTP we present here a complete list of other discovered errors exceeding a 
unit in the last decimal place. 

The original tables (1909) were followed by a flood of corrections, mostly in Arch. Math. 
u. Physik. That perfection was not attained in the later editions is evident. The surviving 
compiler (Emde) is not entirely to blame, as many of the tables have been taken from diverse 
sources (always freely acknowledged); it was obviously impossible to recompute every 
such table. 

Table I isolates errors that occur in the 1909 edition only. Some more known, but ob- 
vious or unimportant, errors have been omitted. In both lists, with few exceptions (the 
principal being in the table of powers in the 1933 edition), errors of a unit in the last decimal 
are not given, although often statistics of their frequency will be found in the notes. The 
practical view-point is that the labour of making such end-figure alterations would not be 
justified; nor is it fair to classify them as errors if the author has not claimed absolute 
accuracy in the last decimal. 

One of my most treasured possessions is the late J. R. Airey’s personal copy of the 1909 
edition of these tables. It is covered with beautifully written notes, references and correc- 
tions, and has 120 pages of related formulae, etc. pasted in—all neatly transcribed, with 
their sources if they are taken from published works. Whenever Airey visited me, he always 
came hugging this book mysteriously! It is now shared with J. C. P. Miller, on whom 
Airey’s mantle has fallen. 

In the date column, . . . denotes that the quantity to be corrected is not given; C de- 
notes that it is correct. Page numbers in square brackets in the 1938 column are a reminder 
that these pages are at the end of the reprint (1943) by Dover Publications. 

In the column “authority,” A = Airey (marked in his personal copy, with no source 
ascribed), C = Comrie, M = Miller, A.M., followed by a number, denotes Archiv fir 
Mathematik und Physik, dritte Reihe, and the volume number. The complete references 
here are: 


v. year Pages authors 

15 1909 372-4 Bruns, Egloff, Emde, Greenhill, Kritaniger, Witt. 
16 1910 285-6 Baruch, Debye, Emde, Goldscheider, Naetsch, Russell. 
17 1910 286 Dinnik, Goldscheider, Jahnke. 

18 1911 204 Dinnik. 

18 1911 374 Emde, Jahnke, Schaefer. 

20 1912 95 Escott. 

21 1913 191 Emde, Escott, Jahnke, Nagaoka. 

22 1913 95 Wagner. 

22 1914 285 Airey. 

22 1914 376 ? 

24 1915 94 Kalahne. 
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TaBLE I-—Errors in 1909 edition only 


for 
xtanx x =0.2 -04045 
x= 1.4 8.2627 
= 1.5 21.1523 
02 1.0112 
x =0.5 1.0927 
x = 1.2 2.1434 
i4 4.2157 
cosx coshx = 1.x = Ois 
also a root 
for x = 6.0 403.00 
Penultimate line —Ar Sin ix 
Line—10 sin 
Last line, middle integral Tes 
1 1 
Tez 
Last line of Section 1 Ei(ix) Ci(x) 
Line 4 — 47 cos 2x 
Ci(x) for x = 40(5)75, 
120(10)170, insert an 
additional cipher after 
the decimal point 
S(1.6) .6383 
C(1.8) .3363 
S(2.3) .5525 
‘ 1 
Line 2 
Line 7 2 
Line 3 of Section 5 See note 
Line 8, denominator x—4%4 
Line 9, denominator 
Line4, replace last term by 
Legend of Fig. 8 Gamma- 
funktion 
Diagram 
Diagram 
Characteristic of logs in 
¢@ = 0°, 5°, 10° 1 
r=48 ¢ = 10° 0 
Line 2, last term in numer- 
ator dv 
Line—3 Amp x 
log g for a = 64° 5’ 4087 
log g for a = 69° 55’ .1195 
Upright side. Denomi- dv 
nators dv, 
Legend of Fig. 22 1.53995 
Beginning of line 7 e 
Line 1 of Section 3; 
Niven’s initials M 


read 
-04054 
8.1170 
21.1521 


1.0136 


1.0926 
2.1435 
4.1413 


Ei(ix) = Ci(x) 
4x(1—cos 2x) 


.6389 
T 
2v2 
x+4 


(3) 


Fakultat 


authority 193 


A.M., 16 
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M 
M 
A.M., 18 
403.43 A.M., 20 
—iArSinix A.M., 16 
Sin 2x A.M., 16 
Tg x A.M., 22 
cm Ctg A.M., 22 
C 22,23 
M 
(12) 
| A.M., 16 
A.M., 16 
A(1) 
A.M., 16 
A.M., 16 
A.M., 16 
A 
36 A.M., 15 
F T A.M., 15 
ee: 45 I M 
i: dnv A.M., 16 
Amp u A.M., 16 
.0087 A.M., 21 
.1159 

c 71 A.M., 16 

iim, 1.52995 A 

dK 

W.D. 


1933-8 


aan 


1909 


87 


111 


125 
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P;(.03) 
P,(1°) 
P;(3°) 
P,(7°) 
P,(7°) 
P;(13°) 
P;(18°) 
P.(19°) 
P,(19°) 
P;(21°) 
P.(25°) 
P;(25°) 
P,(40°) 
P;(40°) 
P,(47°) 
P.(47°) 
P7(47°) 
P,(57°) 
P,(71°) 
P,(71°) 
P,(73°) 
P;(73°) 
P,(74°) 
P,(74°) 
P,(76°) 
P,(76°) 
P;(76°) 
P,(77°) 
P,(77°) 
P;(88°) 


Second formula for Yi(x) J:(x) 


TaBLe I—(Continued) 


for read 
.3000 0300 
.9955 
9617 -9620 
8476 8492 
7986 
3940 .3980 
0289 .0248 
1347 
0443 .0433 
1662 -1664 
-2053 .2040 
3236 
1003 -1006 
4252 4227 
0645 .0665 
2054 -2028 
2949 .2947 
1786 1791 
.1808 
1144 -1136 
2347 .2352 
0795 .0788 
2559 .2563 
4112 -4122 
0076 0070 
.2848 -2850 
0284 .0290 
.2919 
.0735 0755 


lim Ni(+ix) See note (3) 
Expression for large roots 

of Jo 50661 -050661 
Expression for large roots 

of Ji 15399 .015399 
Line — 3 185 191 
Coefficient of cotg*r 

in A2(r) 77 
Coefficient of cotg*r 576 576 

in B2(r) 

1 1 
Series for Jp(x) 4 
Smallest root of J_1/3(x)=0 1.88 1.8663: 
D(10) 078 .068 
Jo(0.62) -9052 -9062 
Jo(5.90) 1120 -1220 
— J,(10.30) + 
Yo' (5.2) .0816 .0846 
Yo’ (8.1) -1806 .1804 
Yo'(8.5) -0091 -0094 
Yo' (8.6) .0336 .0333 


authority 
A.M., 15 
C(2) 


C(2) 


C(2) 


(In x—1) A.M., 16 


A.M., 16 


A.M., 16 
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— 
85 
C 86 
Cc 
96 
Cc 99 
... 100 A 
A.M., 17 
A.M., 17 
A(4) 
A 
| | A.M., 22 
115 A.M., 22 
118 A.M., 16 
A(S) 


1933-8 1909 


127 


128 


C 130 
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TaBLe I—(Continued) 


Yo'(8.7) 

Y'(8.8) 

Yo’ (8.9) 

Yo’ (9.0) 

Yo (9.5) 

Yo'(9.5) 

Yo' (9.6) 

Yo (9.7) 

Yo’ (9.7) 

Yo’ (9.8) 

Yo (9.9) 

Yo'(10.2) 

K,(0.20) 

K,(0.65) 

K;(2.4) 

Ki(8.1) 

K,(8.5) 

Ki(8.6) 

Ki(8.7) 

Ki(8.8) 

K,i(8.9) 

Ki(9.0) 

K,i(9.5) 

K,i(9.6) 

Ko(9.7) 

K,i(9.7) 

Ko(9.9) 

K,(10.2) 

— N,(0.1) 

— N,(0.2) 
No(3.9) 

— N,(5.2) 

— Ni(8.8) 
No(9.9) 

— N,(10.2) 

x = 5.1 to 6.0. Imaginary 
part 

x = 3.1 

x = 3.3 

x = 6.0. Real part of Ho 

x = 2.4, Real part 

Table of ber x, etc. 

J7(13) 

Line—4 

Right-hand side of fifth 
integral 


Last differential equation 


c) Second equation 


Reference to W.D. Niven 
Reference to P. A. Hansen 


1261 
3923 
7.0317 
3.3235 
.0237 
0773 
0541 
.2507 


1.8918 
0.3556 
+ 


See note (8) 


XVI 


See note (9) 
x? 
77 


170 


5.2211 
1.8485 
-2091 


0.8919 


A(6) 


A(6) 


164 
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Sor read authority 901 
0751 0753 9 
1157 1160 10 
.1548 
1921 .1923 
.2463 .2465 
.3378 34 
3574 3571 34 
1755 1752 
3729 3727 
3846 
.0982 .0984 
3931 .3922 
= 5.2209 3s 
1.8483 39 
.2093 74 
0408 0411 
.0020 0017 94 
0438 .0440 103 
0851 .0854 105 
1251 1255 
.1637 .1639 
.3194 .3192 
13412 ‘3410 
.2012 .2009 
3594 
1263 
Bee Cc 129 6.4590 M(7) 
3.3238 M 
0234 M 
.0792 A 
0544 M(7) 
M 
(2502 A 
Cc 137 
+ A.M., 17 
1.3556 A.M., 18 
C A.M., 16 Cc 
+ A.M., 18 
A.M. 
Cc 153 A 
164 | XXVI A 
A.M., 16 = 
A.M., 16 
78 A.M., 20 
= 169 A.M., 15 


MATHEMATICAL TABLES—ERRATA 395 


} TABLE II 
1938 1933 1909 for read authority 
9 86 22 Si(65) 1.5775 1.5792 M(10) 
& 0,067 353011 0,067 352301 A. Fletcher (11) 
571 571 
~ 2488 320x Fletcher 
x 
34 106 ... #,(5) M 
34 109 26 C(0.1) .0999 -1000 (12) 
C(0.8) .7230 
C(1.1) .7638 
C(4.5) .5258 -5260 
C(5.1) 4987 4998 
S(7.3) .5199 .5189 
35 108 25 S(1.5) A153 A155 M 
39 113 ... Heading x 01 A M 
70 142 62 E(15°, 54°) 9345 .9346 M(13) 
74 146 ... Coefficient of 27 1701 1707 A(14) 
94 160 ... Infirst bordered table +4 = Todd (letter) 
103 169 74 ouforr = 123 1.0441 1.0438 M 
105 171 50 Last integral formula for 
Weierstrassian functions. 
After this integral add: 
where v is defined by au au 
e(v) = 8/y + 
108 174 ... Diagram is in error (15) 
117 183 82 Expression for dP;/dé sin 36 sin* @ A.M., 16 (16) 
124 190 88 dP,/d0,n = 3,0 = 31° 2.0654 2.0656 NYMTP (17) 
n=4,@= 5° 0.8570 0.8567 
n= 4,0 = 64° 1.6306 1.6300 
n= 5,0 = 54° 2.0173 2.0178 
n= 5,0 = 64° 1.5420 1.5418 
n = 5,0 = 67° 1.1187 1.1183 
2 = 5,0 = 83° 1.4831 1.4827 
2 = 6,6 = 23° 3.0917 3.0902 
n = 6,0 = 48° 2.3715 2.3713 
n = 6,0 = 57° 1.1936 1.1934 
2 = 6,0 = 83° 1.4487 1.4484 
nm = 7,0 = 70° 1.984 1.934 A 
2233 ... Jae .3250 .3255 M(18) 
Jsi4(4.8) 3518 
J—114(4.8) .1149 
J_si4(4.8) -1803 .1791 
J~14(5.4) .0942 .0940 
J—s14(5.4) .3085 .3082 
164 223 ... Jin(74) .2925 .2923 NYMTP(19) 
.2873 
.2449 
Jsi4(6.8) .0447 
Jsi4(7.8) 


Jaia(8.0) 2749 


\ 


1933 1909 


239 147 


285 


313 


316 
77 
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TaBLe II—(Continued) 


J—114(.2) 

J_314(3.2) 

p=2,n=9 

p=3,n=3 

p=3,n=9 

p=4,n=1 

p=4,n=9 

p=5,n=1 

p=5,n=4 

p=5,n=5 

p=5,2=9 

x = 6.4 

0 and 1. Delete all 
bracketed digits 

(0.20) 

o(10.35) 

Mof13.25) 

$0(4.41) 

$o(8.17) 

Io(6.4) = Jo(6.4:) 

Tyja(1.4) = 4-13 

T113(3.8) 

T13(4.0) 

T_113(3.0) 

I_13(3.6) 

T_113(3.8) 

T_113(4.0) 

T2i3(2.4) 

I2;3(3.6) 

I2/3(3.8) 

I2i2(4.0) 

T2i3(4.4) 

Taia(S.8) 

I_2/3(2.4) 

I_213(2.8) 

I_212(3.4) 

I_2/3(3.6) 

T_2/3(3.8) 

T_213(4.0) 

I_213(4.4) 

T_2/2(4.8) 

T_213(8.0) 

bo(5.3) 

bo(9.5) 

bo(9.6) 

$rbo/bi for r = 0.5 

(4) Comrie’s initials 

(i), line 3 


for 

1.4310 
.3943 
30.571 
13.017 
32.050 
7.586 
33.512 
8.780 
18.982 
22.220 
34.983 
d = 193 


authority 


A, M(20) 


M 


M 
M(21) 
M 


M(21) 

M 

NYMTP 
NYMTP(22) 


1938 


[6] 


[9] 
[15] 


[56] 


(1) 


(2) 


exp 
tha 
ing 

mer 


(3) 


(4) 
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1.4319 
[2] 
168 30.569 
13.015 [4] 
32.065 
33.537 
8.771 [5] 
18.980 
22.218 
34.989 
er 183 253 ... d = 103 = 
8 206 276 163 
135 235 
273 278 
136 236 
eo 228 280 ... 96.98 96.96 
235 1.4002 1.4000 
7.8727 7.8729 
Ric 9.3473 9.3470 
11.1136 11.1138 
2.5226 2.5627 
4.9710 4.7754 
7.8831 7.8827 
9.3546 9.3548 
as 11.1204 11.1201 
65.55 66.55 
a 2.6650 2.6648 
3.7343 3.7345 
7.4423 7.4429 
8.8698 8.8692 
ae 10.5796 10.5799 
15.107 15.098 
53.77 53.71 
2.7065 2.7067 
3.7002 3.7603 
6.2655 6.2657 
= 7.4535 7.4531 
8.8841 8.8773 
10.5867 10.5864 
ea 15.112 15.103 
ee 21.621 21.623 
: 418.01 415.01 A.M., 22 
7.340 7.475 M(23) 
108.02 108.00 
115.31 115.30 
266 | 006 008 M(24) 
300 L.T. LJ. 
.0350 .0400 c = 
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TaBLe II—(Continued) 


1938 1933 1909 for read authority 
(1), line 2 Mittlin Mifflin 
[2] 2... x*forx = 0.83 4747 4746 M(25) 
x® for x = 0.88 4640 4644 
[4] 4 ... x*forx = 1.18 2289 2288 
Powers of x = 1.68 <x® 2222 2248 
x? 3732 3777 
x® 6270 6346 
[5] 10534 10661 
x0 17697 17910 
xu 2973 3009 
4995 5055 
x13 8391 8492 
xia 14097 14267 
xis 2368 2397 
x! for x = 1.40 2892 2893 
{6} 6... fors = 2.05 3620 3621 M(25) 
x? for x = 2.35 3960 3958 
Powers of x = 2.95 x® 6590 6591 
x? 19441 19443 
x® 5735 5736 
[7] 16919 16920 
xl 14723 14725 
x2 4343 4344 
12813 12814 
Powers of x = 2.85 x! 3523 3535 
xu 10041 10076 
x2 2862 2872 
x8 8156 8184 
xia 2324 2333 
xis 6624 6648 
x" forx = 3.10 2540 2541 
x'4 forx = 3.30 18164 18163 
[9] 9 ... x°-forx = 0.1 0.251 0.2512 M 
{15} 15 ... See note (26) 
dx 1 
[56] 56 13 —1 Tex Ctg x A.M., 22 
NotTEs 


(1) Airey has corrected his copy from 


Z (1 — to (1 — + (1 — 7) 
r=2 r=1 

(2) In this table many end-figure changes have been made in the 1933 and 1938 editions, 
expecially in Ps and P;; only those greater than one unit are listed here. It will be noticed 


that, in every case where P¢ is in error, P7 is in error by 2 P, times the error in P¢, thus show- 


ing that P; was found by recurrence, and was not checked by differencing. For further com- 
ments on the original table, see MTAC, p. 137. 
(3) This should read, according to Goldscheider, 
mx” 2 a y¥/ 


(4) See MTAC, p. 196 where the value of the root to 15D is also given. 


@ 
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(5) There are also 44 alterations of a unit in the last decimal, taken from Airey’s table in 
B.A.A\S., Report, 1914. See also Note (6). 

(6) Marked in Airey’s personal copy, together with over 120 alterations of a unit in the last 
figure. See also A.M., v. 18, p. 374 where Jahnke notes these end-figure corrections, based 
on the B.A.A.S., Report, 1911 (computed by Airey). Note the obvious connection between 
the erroneous arguments here and in the values covered by Note (5), which is explained by 


(J.&E.) = 5 (Watson) + In 


K, (J.&E.) = — 3 Y, (Watson) = G, (Airey) 


Watson's Y, is called N, in Jahnke & Emde. 
(7) See the argument 10.2 covered by Notes (5) and (6). Miller has compared this table 
with Watson’s Yo and ¥;. Besides the corrections given, he finds 67 errors of a unit, and 6 of 
two units. Kalahne also points out the error in — N;(0.1) in A.M., v. 24, p. 94. 
(8) See A.M., v. 16 and 21. There are many errors, but it is not worth correcting this table, 
since correct values may be obtained from the later editions (up tox = 10) or from Dwight’s 
Mathematical Tables (up to x = 20). For a thorough discussion of the errors see MTAC, 
p. 297f. 
(9) A further term is needed, namely 

_ (ax) 

(10) Miller has compared these integrals with the NYMTP tables up tox = 100. There are 
also 5 end-figure errors in the 1933 and 1938 versions of this table. 
(11) This error occurs also in Schlémilch, Comp. d. héh. Analysis, v. 2, 1879, p. 258, and in 
Laska, Sammlung, 1888-94, p. 272. A. FLETCHER. 
(12) There have been “errors within errors” here. A first list was given by Lash Miller and 
Gordon in J. Phys. Chem., v. 35, 1931, p. 2872. A reference to this paper, and extracts, were 
given by J. H. Awbery in Phys. So. London, Proc., v. 51, 1939, p. 199. Lash Miller and 
Gordon failed to note the error at S(2.3) and gave C instead of S for the corrections at 1.6 
and 7.3. The first three errors were found independently by Rankin (letter from L. Rosen- 
head, Feb. 1943). In a 4-figure table prepared by Comrie, C(1.8) appears as .3336, and 
S(2.3) as .5532. This table shows that there are no other errors exceeding one unit up to 
u = 3.3. This note has reference to T. I, p. 26 and T. II, p. 34, 109, 26. 
(13) Traced by means of errors in Legendre’s table as given by Heuman (1941). 
(14) Also reported by John Todd by letter. Airey has inserted coefficients up to «** by hand 
in his copy; this one was first written (p. 47, 1909) as 1701 and later changed to 1707. 
(15) G. I. Taylor has pointed out that this diagram is completely wrong near x = — 1. 
Also it should be symmetrical about » = — }, not » = 0, as seems suggested. 
(16) It is surprising that this published correction (which I have verified) was not incor- 
porated in later editions. 
(17) The NYMTP found also 41 errors of a unit in the last decimal. All errors of this table 
are corrected in the sixth printing of Jahnke & Emde. 
(18) Miller has pointed out that the table in 1933 is taken from Dinnik, A.M.., v. 21, 1913, 
p. 324-6, whereas that in 1938 has been corrected to agree with Karas, Z. angew. Math. u. 
Mech., v. 16, 1936, p. 248, who removed some of the major errors in Dinnik and subtabu- 
lated from interval 0.2 to interval 0.1. There are also 5 cases of a change of a unit in the 
last decimal. See MTAC, p. 365-366. 
(19) The NYMTP found also 22 errors of a unit in the last decimal. All errors of this table 
are corrected in the sixth printing of Jahnke & Emde. 
(20) Also 18 errors of a unit in the last decimal. Compare MTAC, p. 160, 282. 
(21) Miller reports that there are also 40 errors of one unit in the last decimal in the table, 
p. 212, 214, 216. This table is theoretically identical with that on p. 218, 220 and 222, 
in which only 8 such end-figure errors have been found, none of which corresponds to any 
of the 40 just mentioned. 
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(22) The NYMTP found also 37 errors of a unit in the last decimal. All errors of this table 
are corrected in the sixth printing of Jahnke & Emde. 

(23) Miller has also found 9 errors of a unit in the column bo and 4 in b;. See also Note (24). 
(24) Miller has also found 17 errors of a unit in the last decimal in the column $rb0/b; up 
tor = 5. See also Note (23). 

(25) This table of powers (for powers greater than 1) has been completely checked by Miller. 
All errors are given here. The errors in the powers of 1.68 are due to taking x* = 1.6625, and 
those in 2.85 to taking x" = 2.84x*. They could have been detected by differencing. 

(26) Note that d is a mean difference, so that the values of the tabulated functions at 1.00 
cannot be accurately obtained; for the upper line of each table this value is 1.0987, and for 


the lower line 0.4551. 


UNPUBLISHED MATHEMATICAL TABLES 
Unpublished mathematical tables are listed in RMT 186 (Lehmer); QR 14. 


32[A].—Wii1aM Pitt DurFEE (1855-1941), Factor Table of the Sixteenth 
Million. Unique ms. calculated during the years 1923-1929, and the 
property of the American Mathematical Society in New York City, 
since December 1935. 


The following description of the ms. was published in Scripta Math., v. 4, 1936, p. 101: 

“The table comprises 500 separate sheets, 8} X 14 inches, each accounting for 2000 
numbers, but as the multiples of 2, 3, 5, and 7 are omitted, the actual entries on each sheet 
number about 416. The entries are in long-hand, in black ink, except that those numbers 
whose lowest prime factor is 11 have been interpolated in red. They are arranged in parallel 
columns, three centuries to a column, the last four digits only of each number being written; 
and opposite each its lowest prime factor. If the number is prime a bar is drawn across the 
corresponding space in the column of prime factors. The arrangement is thus closely similar 
to that in the published tables covering the first nine millions (Burckhardt, Dase and 
Rosenberg, Glaisher). 

“At the foot of each column the number of primes for each century is noted, and the 
total number of primes on preceding sheets, the number on that sheet, and the total. In 
the lower right-hand margin there are listed the number of entries on preceding sheets, 
on that sheet, and the total, with a similar notation in red ink for the entries whose least 
prime factor is 11.” 

The computations, which were made by the stencil method, have not been generally 
checked except by the author of the table. However, D. H. L. discovered one error where 
15485303 is entered as a prime, when 109 is a factor. The total number of primes in this 
million is 60,465. 


33[A].—U. S. Dept. or ComMERCE, WEATHER BurREAU, Table of (log e)/x, 
computed by, and in possession of, the Bureau. 
The Table of (log e)/x is for x = [100(1)999; 6S]. 
L. P. HARRISON 
Weather Bureau 


MECHANICAL AIDS TO COMPUTATION 


14[Z].—G. W. Kine, ‘“‘Punched-card tables of the exponential functions.” 

Rev. Sci. Instruments, v. 15, 1944, p. 349-350. 

G. B. THomas & G. W. Kine, “Preparation of Punched-card tables of 
logarithms,” ibid., p. 350. 

With the rapid increase in recent years of the use of punched-card machines in com- 
putations of all kinds, the versatility of the machines has been adapted to many uses. 
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One obvious method of increasing the number of possible applications has been to prepare 


cards punched with the values of numerical functions which can be combined with other 
values by multiplication or addition. Of this type are exponential and logarithmic functions. 

The first note indicates the usefulness of punched-card tables of the exponential func- 
tion. From the relationship, exp (a + b) = (exp a)(exp 5), one sees that this function is 
readily adapted to calculation by punched-card machines, any desired degree of accuracy 
being easily attained. Thus, the master set of cards might contain values of exp + n over 
the range m = 0.00(0.01)12.00. A second set of 100 cards for n = 0.0000(0.0001) (0.0099) 
would then multiply the original range a hundredfold. 

One of the most important problems in the use of punched cards is the verification of 
the numbers punched on them. For this purpose the author suggests the use of the formula 


Partial sums of the cards are compared with the values of this function computed for 
suitably chosen values of Ak. 

The second note suggests the use of punched-card tables of log x. The authors suggest 
that the verification of the cards be made by comparing partial sums of the cards with 
computed values of log m! from the formula: 


DY log & = log n!. 
k=1 


Tables such as those of J. Boccarpt: Tables logarithmiques des factorielles jusqu’a 10,0001, 
Cavaillon, 1932, or F. J. DUARTE: Nouvelles tables logarithmiques @ 36 decimales, Paris, 1933, 
are available for this purpose. 

T. BD. 


NOTES 


29. EaRLy DeciMAL DIVISION OF THE SEXAGESIMAL DEGREE.—In 
MTAC, p. 33, 100 (corrigenda), 129-130, it has been already noted that 
decimal division of the degree was used by RuFF! in a Latin codex of about 
1450, as well as by BricGs, in tables prepared before 1633, the idea having 
been suggested to him by Viéte in a work published in 1600. The decimal 
division of the degree was also advocated by Simon StTEvIN (1548-1620) 
of Bruges in 1585, in his Flemish work, 1. De Thiende, published at Leiden. 
A facsimile of this 36-page booklet is given in H. Bosmans, La ‘Thiende’ 
de Simon Stevin, Antwerp and The Hague, 1924; a French edition of De 
Thiende, La Disme, occupied p. 132-160 of Stevin’s 2. La Pratiquve D’ Arith- 
metiqve, also published at Leiden in 1585; there are copies of this work in 
the Plimpton Library of Columbia University, the Royal Library of Bel- 
gium, University of Liége, etc. A facsimile of this edition was published 
by G. SaRTON in “The first explanation of decimal fractions and measures 
(1585). Together with a history of the decimal idea and a facsimile (no. 
XVII) of Stevin’s Disme,”’ Isis, v. 23, 1935, p. 230-244. A French edition, 
ed. by A. Girard, appeared also in a new edition of 3. L’ Arithmetique, Leiden, 
1625, p. 823-849, and in 4. S. Stevin, Les Oeuvres Mathematiques, Leiden, 
1634, p. 206-213. There were other Flemish editions, 5—6, in 1626 (in the 
Brown University Library), and in 1630. 

This was the first work to set forth the theory of decimal fractions. It 
was translated into English by 7. RopERT Norton, “engineer and gunner” 
(d. 1635, Dict. Nat. Biog.), son of the poet THomas Norton (1532-1584), 
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under the title Disme: The Art of Tenths, or, Decimail Arithmetike,Teaching 
how to performe all Computations whatsoeuer, by whole Numbers without 
Fractions, by the foure Principles of Common Arithmeticke . . . Inuented 
by the excellent Mathematician Simon Steuin, London, 1608; there are copies 
of this work in the British Museum, the Cambridge University Library, and 
in the Harvard University Library (a copy formerly in the W. A. White 
Library). According to the Cambridge Bibliography of English Literature, 
there were other editions, 8(?) in 1614 and, by HENry LyTE, in 1619. The 
last statement is certainly incorrect in the case of Lyte’s independent work 
(with a similar title), of which copies are in the Huntington Library, Pasa- 
dena, and the British Museum. See further Sarton, Lc., p. 158f. Norton 
supplied tables of interest and measurement, and instructions in decimal 
arithmetic to Robert Recorde’s The Grovnd of Artes, London, 1623. See 
also H. Bosmans, ‘‘La ‘Thiende’ de Simon Stevin,” Rev. d. Questions Sci., 
s. 3, v. 27, 1920, p. 109-139. 

The first of two paragraphs in ‘‘Article V, Des Compvtations Astro- 
nomiqves’”’ of the ‘““Appendice” of La Disme (1585), with its illustration of 
Stevin’s decimal notation (by ‘“‘commencemens” he referred to the integral 
part of the decimal) is as follows: 

Aians les anciens Astronomes parti le circle en 360 degrez, ils voioient 
que les computations Astronomiques d’icelles, auec leurs partitions, estoient 
trop labourieuses, pourtant ils ont parti chasque degré en certaines parties, 
& les mesmes autrefois en autant, &c. a fin de pouuoir par ainsi tousiours 
operer par nombres entiers, en choissisans la soixantiesme progression, parce 
que 60 est nombre mesurable par plusieurs mesures entieres, 4 scauoir 
1, 2, 3, 4, 5, 6, 10, 12, 15, 20, 30, mais si l’on peut croire l’experience (ce 
que nous disons par toute reuerence de la venerable antiquité & esmeu 
auec |’vtilité commune) certes la soixantiesme progression n’estoit pas la 
plus commode, au moins entre celles qui consistoient potentiellement en la 
nature, ains la dixiesme qui est telle: Nous nommons les 360 degrez aussi 
commencemens, les denotans ainsi 360(0), & chascun degré ou 1(0) se diuisera 
en 10 parties egales, desquelles chascune sera 1(1), puis chasque 1(1) en 
10(2), & ainsi des autres, comme le semblable est faict par plusieurs fois ci 
deuant. 

[In 4 places near the end of this quotation we have, for typographical 
reasons, printed (0), (1), (2), instead of the original notation of the figures 
0, 1, 2 surrounded by circles.] 

R. C. A. 


30. First PuBLIsHED CoMPOUND INTEREST TABLES.—These were 
1. Simon Stevin’s Tafelen van Interest, Midtsgaders De Constructie der seluer, 
[Tables of interest together with their construction], Antwerp, 1582, 92 p. 
A French translation from the Dutch was published in 2. Stevin’s La 
Pratiqve D’ Arithmetiqve, Leiden, 1585, and an edition 3. by Girard in 1625 
was reprinted in 4. Stevin’s Oeuvres Mathematiques, Leyden, 1634, p. 185— 
206. A facsimile of the original work is given in 5. C. M. WALLER ZEPER, 
De oudste intresttafels in Itahié, Frankrijk en Nederland met een herdruk van 
Stevins “‘Tafelen van Interest,’’ Diss. Leiden, Amsterdam, 1937. 

Stevin tells us in his preface (dated Leiden July 16, 1582) that the 
inventor of these tables was JEAN TRENCHANT, who gave specimens of them 
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in his L’Arithmetiquve de Ian Trenchant departie en troys livres, Ensemble un 
petit discours des changes, auec L’art de calculer aux Getons, Lyons, 1558. 
See G. Sarton, “JEAN TRENCHANT, French mathematician of the second 
half of the sixteenth century,” Isis, v. 21, 1934, p. 207-209; and H. Bos- 
MANS, Soc. Sci. de Bruxelles, Annales, v. 33, 1909, part 1, p. 184-192. 

Since the description by A. DeMorgan of Stevin’s Tables (Supplement 
to the Penny Cyclopaedia, v. 2, London, 1846, p. 604) contains an interesting 
suggestion we shall quote this: “it is pretty certain that these tables of 
compound interest suggested decimal fractions, the account of which speedily 
follows them. They are constructed as follows:—Ten millions being taken 
as the base (or root, as Stevinus calls it), and a rate, say five per cent., being 
chosen, the present value of ten millions due at the end of 1, 2, &c., up to 
30 years, are put in a column, to the nearest integer. By their sides are 
sums of their values, which give the present values of the several annuities 
of ten million, as follows: 


Table d’Interest de 5 pour 100. 


1 9523810 9523810 
2 9070295 18594105 
3 8638376 27232481 
4 8227025 35459506 

30 +«-2313774 


The rates are from 1 to 16 per cent., and also for 1 in 15, 1 in 16, &c., to1 
in 22;? or, as the French say, denier quinze, denier seize, &c. At the end is 
a direction to dispense, when convenient, with some of the last figures.” 

“There is thus a virtual use of decimal fractions preceding the formal 
one. The same thing happens in the tables of RicHarD Wirt, [Arithmeticall 
questions, touching the Buying or Exchange of Annuities . . . , London, 
1613] . . . which we believe to be the first English tables of compound 
interest, and the first English work (except the translation of the Disme 
of Stevinus) in which decimals were used.” 

R. C. A. 


1 Apparently De Morgan is here quoting from Stevin’s Oeuvres. It should be 153724494 
as in the 1582 and 1585 editions.—EpiTor. 


2 The table for 1 in 20 is not given since it would be the same as the 5% table.—Eprtor. 


31. First PuBLISHED MorTALITy TABLE.—At a time when the city of 
Breslau is so much in the war news it is appropriate here to recall that 
mortality records of this city were the basis of a paper, written by EpMOND 
HALLEY (1656-1742) over 250 years ago, and published in R. So. London, 
Phil. Trans., 1693, no. 196, p. 596-610, and no. 198, p. 654-656. It is en- 
titled ‘‘An estimate of the degrees of the mortality of mankind, drawn from 
curious tables of the births and funerals at the city of Breslaw; with an 
attempt to ascertain the price of annuities upon lives.’”’ Halley reprinted 
this in his Miscellanea Curiosa, London, v. 1, 1705; second ed. 1708; third ed., 
1726, p. 280—301. This was also reprinted in Institute of Actuaries, J., v. 18, 
1874, p. 251f and 262f. Here we have the beginnings of statistical science and 
the first published mortality table, in the form since followed, to show the 
annual movement of a population, the problem of survivorship, the average 
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duration of one or more lives, and the money values depending thereon. 
See C. WALFoRD, Insurance Cyclopaedia, London, v. 5, 1878, p. 616-618; 
v. 1, 1871, “Breslau table of mortality,” p. 360-363; H. Béckh, “Halley 
als Statistiker. Zur Feier des zweihundertjahrigen Bestehens von Halleys 
Sterblichkeitstafel,’’ L’Institut International de Statistique, Bull., v. 7, 
1893, p. 1-24; A. Loewy, Versicherungslexikon, ed. by A. Manes, third 
ed., Berlin 1930, cols. 742-744, 998-1000; E. J. Farren, “On indirect 
methods of acquiring knowledge.—The method of history.—The first table 
of mortality,” Assurance Mag., v. 1, 1851, p. 40f. 

Perhaps the inclusion of the following item is warranted, in not being 
too far removed from the main topic of the note. In Math. Gazette, v. 14, 
Dec. 1929, p. 574f, I published a Note on “A letter of DeMoivre and a 
theorem of Halley.”” The Note consisted mainly of a letter from De Moivre 
(1667-1754) to Halley which I copied from The Mathematical Magazine 
and Philosophical Repository, 1761 (see Math. Gazetie, v. 14, p. 387, and 
v. 15, p. 215f) where it is stated that the letter, which is undated, was never 
before published. The theorem gives the rate of interest for an annuity, 
in a compact formula, 1 + 6 — Vé? — 2by (for explanation of the notation, 
see the source). I asked, ‘‘Can any reader state where Halley’s theorem is 
to be found or otherwise assist in determining an approximate date for the 
letter?’’ A long and very interesting reply from F. P. White, of St. John’s 
College, Cambridge, appeared in the Math. Gazette, v. 15, Oct. 1930, p. 
213f. The formula was first given by Halley in his chapter, “Of compound 
interest and annuities,” in the first edition (1705) of Sherwin’s Mathematical 
Tables. Hence De Moivre’s letter, which was found among Halley’s papers 
after his death in 1742, may have been written before 1706. 

Halley was in 1703 the successor of Wallis as Savilian professor of 
geometry at Oxford, and in 1721 followed Flamsteed as astronomer royal. 
Due to his efforts the publication in 1687 of Newton’s Principia was 
achieved. Among his contributions to astronomy was the discovery in 1682 
of the comet, with a period of 76 years, last seen in 1910. He was also the 
translator from the Arabic and restorer of works of Apollonius of Perga 
(Oxford, 1706-1710). 


1It was reprinted in Assurance Mag. and J., v. 9, 1861, p. 259-269. 


32. Guiwe (MTAC no. 7), SuprL. 1.—In F. Bercer, “Uber den Tem- 
peraturverlauf in einem Zylinder von endlicher Lange beim Abkiihlen und 
Erwarmen,” Z. f. angew. Math. u. Mech., v. 11, 1931, p. 53, there are 8 
values (multiplied by 80), to 4—5D, of roots of the equation xJ;(x) = .53Jo(x). 
On p. 50 are graphical solutions of the first five roots. In GERALD PICKETT, 
“The effect of Biot’s modulus on transient thermal stresses in concrete 
cylinders,’”’ Am. So. Testing Materials, Trans., v. 39, 1939, p. 916, a table 
is given of the first six roots of 


xJi(x) = BJo(x), 


B = .1,.5, 1, 3, 10, 100, 1000, 10 000. The first five roots are to 3D and the 
sixth to 2D. 


In H. A. Etnstetn, “Der Geschiebbetrieb als Wahrscheinlichkeits- 
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problem,” diss. Eidgen. techn. Hochschule, Ziirich, Versuchsanstalt fiir 
Wasserbau, Mitteilung, 1937, p. 93, there are graphs of the functions e~*Jo(x) 
and e~*J,(x) on a logarithmic scale in which x ranges from .1 to about 100. 
T. II gives 3D for e“Io(u), u = 2(x7)!, and 4D for exp. [— (x! — T*)*]Je-“Jo(u), 
T = 0.895, 1.689, 3.14, 6.05, 11.48, 18.57; Vx has various values from 0 to 
6.18, a different set for each value of T. 

In I. E. Garrick, “On some reciprocal relations in the theory of non- 
stationary flows,” U. S. Nat. adv. Comm. f. Aeronautics (NACA), Techn. 
Reports, no. 629, 1938, p. 348, for Theodorsen’s function, 


C(k) = + (k)] = + 1G(k), 
there are graphs of 


F(k) = (Ji(Ji + Yo) + Vi — + Yo)? + (Yi — 
G(R) [VY Yo + Ji + Yo)? + (Yi 


0 < k < 2, and also of ¥ (F? + G) and tan (G/F), 0 < k < 1. A table 
with 4D for k = .1(.1).6(.2)1, 2, 4, 6, 10 is given by T. THEODORSEN in the 
paper ‘‘General theory of aerodynamic instability and the mechanism of 
flutter,” NACA, Techn. Reports, no. 496, 1934, p. 426, Jo(k), Ji(k), Yo(k), 
Y,(k), F(R), and with 34D for F, and G, for the same values of k and also 
k = 0, .025, .05, ~. On p. 418 there are graphs of F and G against 1/k 
for 0 < 1/k < 40. On p. 350 Garrick gives a table with 4D of Wagner's 
function k,(x) for x = 0(.5)1(1)10, 20, ©, and also a graph (p. 347), where 


ki(x) = (2/7) sin (xt) F(t)dt/t = 1 + (2/m) cos (xt)G(t)dt/t. These 


values of k(x) are taken from a paper by H. Wacner, ‘Uber die Ent- 
stehung des dynamischen Auftriebes von Tragfliigeln,” Z. angew. Math. u. 
Mech., v. 5, 1925, p. 31. 

In Max Jaco, “Influence of nonuniform development of heat upon the 
temperature distribution in electrical coils and similar heat sources of simple 
form,” A.S.M.E., Trans., v. 65, 1943, p. 597, there is a table of the function 


Jo(xy) — Jo(x) x = 0, .1, 1, 2, 2.3, 2.4048; 
1 — Jo(x) y = 0(.1)1. 


In Fig. 3 the values x = 0, 1, 2 and 2.4048 are used. 

In C. WEBER, “Zum Zerfall eines Fliissigkeitsstrahles,”’ Z. f. angew. 
Math. u. Mech., v. 11, 1931, p. 150, values are given for fo(x), f:(x), where 
= Ko(x)/Ki(x), and fi(x) = x/[xfo(x) + 1], x = [.1(.1)1(.2)2(.5)4(1)9; 


H. B. 


33. A NEw REsuLT CONCERNING A MERSENNE NuUMBER.—On 2 De- 
cember 1944 I finished the proof that the Mersenne number Mg; is com- 
posite; see MTAC, p. 333. Hence there now remain only four of these 
numbers M, (p = 193, 199, 227, 229) whose character is unknown. The 
method of testing primality was the same as that described in Nat. Acad. 
Sci., Proc., v. 30, Oct. 1944, p. 314-316. 


H. S. UHLER 
12 Hawthorne Ave., Hamden 14, Conn. 


| 
34. 
MuIR 
use of 
first p 
notebc 
is use 
radial 
tes of A. 
form 
a the te 
statec 
From 
wholl 
above 
| 
that | 
June 
202). 
: plate 
first 
for 
for 7 
| k-th 
3D]. 
a is, f 
for | 
cien 


fiir 
(x) 
00. 
u), 
| to 


on- 
hn. 


NOTES 405 


34. ORIGINATORS OF THE TERM RADIAN.—As long ago as 1910 THomas 
Murr pointed out (Nature, v. 83, p. 156) that while the earliest recorded 
use of the term Radian in the New English Dictionary was in 1879, in the 
first part of the new edition of the first volume of William Thomson and 
P. G. Tait’s Treatise on Natural Philosophy, ‘‘my own first use of it was in 
class-teaching in the College Hall at St. Andrews in 1869, and I possess a 
notebook, belonging to one of my students of that year, in which the word 
is used.” He hesitated, however, in a definite choice between the terms 
radial, radian, rad. But he states that as a result of reading a publication 
of A. J. Exits (1814-1890) and exchanging letters with him in 1874 “the 
form radian was definitely adopted by me.”’ Ellis remarks that he had used 
the term ‘‘Radial angle” from his Cambridge undergraduate days, but Muir 
stated that Ellis approved of radian as a contraction of “‘radi-al an-gle.” 
From later correspondence in Nature, v. 83, p. 217, 459-460, it appears that, 
wholly independent of Muir, JAMES THomson (1822-1892), brother of the 
above-mentioned William, proposed the name Radian in July 1871 and 
that he used it in an examination paper at Queen’s College, Belfast, on 
June 5, 1873, published in the college calendar for 1873-74. 

Bibliographic reports on the use before 1869 of the term Radial Angle, 
as equivalent to Radian, are desired. This term is not listed in N.E.D. 

A. 


35. Phil. Mag. Tasies, Supp. 3 (for Suppl. 1-2, see MTAC, p. 201-— 
202).—W._G. BickLey, ‘‘Deflexions and vibrations of a circular elastic 
plate under tension,” s. 7, v. 15, 1933, p. 795. The table gives, to 5S, the 
first two roots of 


xIn41(x) V(x? + V(x? + 
+ 
for n = 1,c = 0, 1, 2, 5, 10, 20; the values of x and x*(x* + c*) are given, and 


for n = 2 the same quantities are given for the first root. This item was 
overlooked in the Guide, MTAC 7. 


H. B. 


36. ZEROS OF THE BESSEL FunctTION J,(x).—If we denote, as usual, the 
k-th positive zero of J,(x) by j,,, then the symmetric function 


= 


is, for each positive integer , a rational function of v. It was first used by 
RAYLEIGH! for the calculation of jo: and j:,: and later by ArrEy* and others 
for many values of j,,1. These functions o2,(v) are also important as coeffi- 
cients of the meromorphic functions 
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This last expansion, in effect, was given as far as nm = 4 by JacosI’ in 1849. 
This is doubtless the first appearance of these functions o. Later writers 
have given lists of these functions as follows: 


author range of » 
RAYLEIGH! (1874) 1(1)5, 8 
Grar & GuBLER‘ (1898) 1(1)5 
NIELSEN® 1904) 1(1)5 
KAPTEYN® (1906) 1(1)6 
ForsytH’ (1920) 1(1)3 
Watson® (1922) 1(1)5, 8 


As a by-product of a recent investigation the first dozen of the functions 
were computed and are given below. They have interesting properties which 
may be discussed elsewhere. We nced only the following explanations here. 
If we use [x] to denote, as usual, the greatest integer <x, and if we define 
the polynomial z,(v) by 


wa(v) = (& — 


k=1 


then the function ¢,(v) defined by 
= 
is a polynomial of degree 
[=]. 
k=2 k 
If we write 


then the coefficients a, for m < 12 are given in the following table. 


b n=1 2 3 4 5 6 7 8 9 
0 1 1 2 11 38 946 4580 202738 3786092 
1 5 14 1026 4324 311387 6425694 
2 362 1316 185430 4434158 
3 42 132 53752 1596148 
4 7640 317136 
5 429 33134 
6 1430 
n= 10 il 12 
h 
0 261868876 1992367192 2381255244240 
1 579783114 4028104212 7315072725560 
2 547167306 3458238276 10093635442028 
3 287834558 1649756012 8275251041478 
4 92481350 479550668 4491836314618 
5 18631334 87264812 1701744523728 
6 2305702 9748732 461790600920 
7 160850 614228 90534103564 
8 4862 16796 12743301316 
9 1257916872 
10 82812812 
11 3271518 
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Thus for = 6 we have 


42% + + 1026» + 946 d 
+ 1)%(v + 2)%(» + 3)*%(» + 4)(» + 5)(v + 6) 


D. H. L. 


? RAYLEIGH, London Math. So. Proc., s. 1, v. 5, 1874, p. 119-124; Scient ye 


‘a: 
-_. Roos 19s. The entry for n=8 is due to A. CAYLEY; see also his Rea 
v. 7, 


2J.R. } el Phil. Mag., s. 6, v. 41, 1921, p. 200-203. 

3C. G. J. Jacost, Asir. Nachrichten, v. 28, ‘1849, cols. 93-94; Gesammelte Werke, Berlin, 
v. y, 1891, p. 173 [for 101 + 32, read 101 + 22). 

. H. & E. Gusier, Einleitung in die Theorie der Bessel’ schen Funktionen, v. 1, 

“The 1898, p. 130. 

5N. NIELSEN, Handbuch der Theorie der Cylinderfunktionen, Leipzig, 1904, p. 360. 

SW. KAPTEYN, Archives Néerlandaises d. Sci. exactes et nat., s. 2, v. 11, 1906, p. 149, 168. 

TALR. ForsytH, Mess. Math., s. 2, v. 50, 1920, p. 135. 
leon ia Watson, A Treatise on the Theory of Bessel Functions, Cambridge, 1922 or 

, p. 502. 
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QUERIES 
13. TABLES OF INTEGRALS.—We are now interested in evaluating inte- 


grals of the following forms: f e~‘dt/t", la e~“dt/#". Are there published 


tables of these functions? MELviIn Mooney 


U. S. Rubber Co., 
Research and Technical Development Dept., 
Passaic, N. J. 


EprtoriaL Note: Among many tables of Lf e~‘dt/t = — Ei(—x) reference may 


be given to NYMTP, Tables of Sine, Cosine and Exponential Integrals, 2v., 1940, for 
x = [0(.0001) 1.9999; 9D}, [0(.001) 10; 9S], [10(.1)15; 14D]. There are useful Bibliographies 


in the volumes. When 2 is a positive integer Lf e~‘dt/t* may be made to depend upon 


Ei(—<x). For the cases m = + 2(—1) — 2 tables were published by W. L. Miller & T. R. 
Rosebrugh, R. So. Canada, Proc. and Trans., series 2, section III, v. 9, 1903, p. 80-101, for 


x = [.1(.001)1(.01)2; 9D]. There are also tables (p. 80-81) of — -f e~‘dt/® + 1/x + Inx, 
and — f e~tdt/t — In x, for x = [0(.001).1; 9D]. In the case of the second integral, when 
n = 0 we have the error function of which the most extensive table is that of A. A. MARKov, 
Table des Valeurs de I'Intégrale {- e~®dt, St. Petersburg, 1888, for x = [0(.001)3(.01)4.8; 
11D] with A*; see MTAC, p. 136. However a more extensive table of the closely related 
function H(x) = - f e~*dt has been published in NYMTP, Tables of Probability Func- 

tions, v. 1, 1941, x = [0(.0001)1(.001)5.6; 15D]. This table can be used to evaluate the 
above integral by means of the relation f™ e~*dt = # {1 — H(x)]. Are there other tables 
of the first function than for — 2 > m > 2, and of the second for nm + 0? 
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14. TaBLes or N*? (Q5, p. 131; QR 8, p. 204; 11, p. 336; 13, p. 375).— 
We have ms. tables, to 10S, as follows for: 


N = 100(1)1000, 1000(10)10 000, 1005(10)1565, and also 
N = [1.0001(.0001) 1.0099; 9D]. 


- 
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On setting (1 + a)? = 1 + (3/2)a + ¢, we found that 
for 0 < a < 0.000 036 514, @ = 0.000 000 000; 
for 0.000 063 246 < a < 0.000 081 650, @ = 0.000 000 002; 
for 0.000 096 610 < a < 0.000 109 544, » = 0.000 000 004. 


IRWIN ROMAN 
U. S. Bureau of Mines, 
Custom House, Baltimore, Md. 


CORRIGENDA 
1. 15, for 220, read 120. 
4, for iHo(—ix), read —iHy®(—ix); for Hi®(—ix), read —Hy®(—ix). 


10, for —.4(.01)+1.49, read —.49(.01)+1.49; for .4(.01)+.99, read .49(.01).99. 
33, for +e*w,(—x), read —e*w,(—x). 


1. 11, for e = 0(12)999, read e = cos 0, = 0(19)992, 


P. 168, 
P. 229, 
P 
| 
| 
; 
. 
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